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ABSTRACT 

Forced vibrations in turbomachinery components can 
cause blades to crack or fail due to high-cycle fatigue. Such 
forced response problems will become more pronounced in 
newer engines with higher pressure ratios and smaller axial gap 
between blade rows. An accurate numerical prediction of the 
unsteady aerodynamics phenomena that cause resonant forced 
vibrations is increasingly important to designers. Validation of 
the computational fluid dynamics (CFD) codes used to model 
the unsteady aerodynamic excitations is necessary before these 
codes can be used with confidence. Recently published 
benchmark data, including unsteady pressures and vibratory 
strains, for a high-pressure turbine stage makes such code 
validation possible. In the present work, a three dimensional, 
unsteady, multi blade-row, Reynolds- Averaged Navier Stokes 
code is applied to a turbine stage that was recently tested in a 
short duration test facility. Two configurations with three 
operating conditions corresponding to modes 2, 3, and 4 
crossings on the Campbell diagram are analyzed. Unsteady 
pressures on the rotor surface are compared with data. 

INTRODUCTION 

Forced vibrations in turbomachinery components can 
cause High Cycle Fatigue (HCF) failures in blades. Modem 
turbomachinery blade designs have higher pressure ratios and 


smaller axial gap between blade rows, which can cause an 
increase in forced response problems. These problems are 
typically not discovered until late in the engine development 
stage. This can lead to increased development time and cost, 
increased maintenance cost, and reduced operational readiness. 
Accurate numerical prediction of the unsteady aerodynamics 
phenomena that cause resonant forced vibrations can lead to an 
improved understanding of the problem and offer potential 
approaches to reduce or eliminate specific forced response 
problems. 

Although the Campbell diagram can identify regions of 
potential problems, the actual magnitude of the aerodynamic 
excitation can only be obtained from experimental 
measurements or unsteady aerodynamic modeling. Without 
knowing the level of excitation, it is impossible to determine 
whether a particular crossing will result in problems. With 
recent advances in CFD modeling techniques and increased 
availability of computational resources, it is now possible to 
model a turbine stage in considerable detail. Such modeling, 
which includes three dimensional and viscous effects, can 
provide details of the unsteady flow phenomena and yield the 
aerodynamic excitation to calculate the blade vibratory 
response. However, validation of such CFD-based models 
requires unsteady aerodynamic data. Recently, such data has 
become available for a full-scale turbine stage operating at 
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engine conditions [1]. This data, acquired in a short duration 
test facility (shock tunnel), includes measurements of the 
unsteady pressure on the blade surface and the vibratory 
response of the blades to excitations from upstream vanes at 
different operating speeds. Further, this data set includes 
measurements at two different axial gaps between the stator 
and the rotor. The test article is an initial design for a high- 
pressure turbine stage, which experienced forced response 
vibrations. An increase in the axial gap was required to 
eliminate the forced vibrations. This benchmark data is suitable 
for validation of unsteady aerodynamics and aeroelastic 
analysis codes. 

In the present work, a 3D, unsteady, multi blade-row, 
Reynolds-Averaged Navier Stokes (RANS) code, TURBO, 
developed by Chen et al. [2-4] is applied to the high pressure 
turbine (HPT) stage for which benchmark experimental data is 
available. Two configurations are analyzed with different axial 
spacing between blade rows. Three operating conditions 
corresponding to modes 2, 3, and 4 resonant crossings on the 
Campbell diagram are analyzed. Two grids are used to evaluate 
the effects of spatial discretization. Unsteady blade surface 
pressures are Fourier transformed to time average and first 
harmonic components and compared with the data at 50% and 
85% span, which were the locations of the pressure transducers 
in the experiment. In addition, plots of the flowfield are 
presented to view how upstream wakes interact with 
downstream rotor. The results presented show the effect of 
axial spacing on the unsteady flowfield on the downstream 
rotor. The paper summarizes the results of the comparison of 
the numerical predictions with data, and concludes with some 
recommendations for future work. 

DESCRIPTION OF EXPERIMENTAL DATA 

The detailed experimental data was published by Kielb et 
al. in Refs. [1,5,6]. The data was acquired in a short duration 
test facility at the Gas Turbine Lab of the Ohio State 
University. The test article was a high pressure turbine in which 
the rotor has experienced excessive vibrations due to pressure 
wakes from upstream stators. The on-blade instrumentation 
included pressure transducers and strain gages. The data was 
acquired at and near the resonant crossings of the second (1 st 
torsion), third (axial) and fourth (2 nd bending) modes with the 
41 engine order excitations from upstream vanes. 
Measurements were obtained for two axial spacings between 
the stator and rotor. The surface pressure measurements were 
presented as plots of time-averaged and two harmonics 
components with axial chord location at 50% and 85% span 
locations. 

DESCRIPTION OF TURBOMACHINERY UNSTEADY 
FLOW CODE 

The Reynolds-averaged Navier Stokes turbomachinery 
code used in the present study is the TURBO code [2-4, 7, 8]. 
This 3D, unsteady, multi blade-row code is a structured grid 


solver based on an implicit finite volume algorithm. Flux 
vector splitting is used to evaluate the flux Jacobians on the left 
hand side of the governing equations and Roe's flux difference 
splitting is used to form a higher-order TVD (Total Variation 
Diminishing) scheme to evaluate the fluxes on the right hand 
side. Newton sub-iterations are used at each time step to 
maintain higher accuracy. Symmetric Gauss-Seidel iterations 
are applied to the discretized equations. Viscous terms are 
modeled explicitly and the two-equation k-epsilon turbulence 
model is used for closure. 

For rotor-stator interaction calculations, various 
approaches are possible to model the unsteady flow field. If 
periodic boundary conditions are used in the blade-to-blade 
direction, it may be necessary to include multiple blade 
passages in each blade row. For example, in the present work, 
the turbine has 41 stators and 78 rotors. Use of periodic 
boundary conditions would require calculations to include 41 
stator passages and 78 rotor passages. Alternately, one could 
modify the number of stators from 41 to 39 and change the 
solidity accordingly so that the calculations can be done with 1 
stator passage and 2 rotor passages. 

The approach adopted in the present work is to model the 
actual geometry without modifications and to use phase-lag 
boundary conditions so that only one passage from each row 
needs to be included in the calculations. The boundary 
conditions are based on the direct-store method [9] and the 
implementation of these in the TURBO code is described in 
Refs. [3,4]. The history of flow variables is stored and 
continually updated at the interface between the two blade 
rows. These stored flow variables are available as interface 
boundary conditions for the rotor and stator blade passage, 
regardless of its relative location. The same technique is 
applied to the circumferentially periodic boundaries of both 
rotor and stator passages. 

DETAILS OF UNSTEADY FLOW COMPUTATIONS 

The computational grids used in the present calculations 
are smoothed H-grids. The stator and rotor passage grids were 
95x43x37 (axial x radial x circumferential) and 93x43x37, 
respectively, with 55 points in the axial direction on each 
surface between leading and trailing edges. Figure la and lb 
show sample grids at the mid-span location in the cascade 
plane. Although multiple passages are shown for clarity, the 
calculation domain includes only one blade passage in each 
blade row. 

The computations using the TURBO code were started 
from a uniform flow and converged to periodicity. 
Computations were performed for two values of axial spacing 
between stator and rotor: 32% and 46% of stator chord, 
referred to as the small and large gap. The rotational speed, 
inlet total conditions, and exit static pressure used in the 
TURBO calculations were the values reported in Ref. [6]. Table 
1 summarizes the conditions for the 6 different runs. The run 
numbers listed correspond to those in Ref. [6]. 
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RESULTS 

Initial results presented in this section seek to demonstrate 
the periodicity of the unsteady solutions, and the effect of the 
time step, and other numerical parameters on the unsteady 
pressure results. These calculations were performed for the Run 
1 conditions. The results presented later in this section show 
the comparison between the computational results and 
measurements at two spanwise locations for Run 1 and 11 
conditions. Finally, select results are presented for Run 4 and 7 
conditions. Note that surface pressures have been normalized 
by the inlet total pressure. 

The periodicity of the flowfield was checked by 
monitoring the unsteady pressure at many locations. In Fig. 2, 
the variation of unsteady blade surface pressure at selected 
locations is plotted against time. The quantity plotted is the 
difference between instantaneous pressure and the time- 
averaged pressure, scaled by the inlet total pressure. The plots 
are for the two spanwise locations at which measurements were 
obtained. The chordwise locations are near the leading edge, 
near the mid-chord, and near the trailing edge on both blade 
surfaces. 

It can be seen that excellent flow periodicity is present in 
the vicinity of the leading edge and mid-chord. Near the 
trailing edge location at 50% span, small differences are seen in 
the surface pressure from one cycle to the next. Near the 
trailing edge location at 85% span, the lack of periodicity is 
quite noticeable. Note that this aperiodicity is restricted to a 
small region near the trailing edge and is not eliminated by 
continuing the unsteady calculations for a longer time. The 
dominant frequency of the unsteadiness in the trailing edge 
region is not the blade passing frequency, but nearly six times 
larger. These high frequency variations in pressure are caused 
by vortex shedding at the rotor trailing edge, which is clearly 
observed in the flow visualization plots presented later (for 
example, Fig. 5). The difference in periodicity between the 
85% and 50% span locations is also caused by differences in 
the vortex shedding patterns as seen in these flow visualization 
plots. 

Effect of Numerical Parameters 

In Fig. 3, the effect of time step size on the unsteady 
pressure results is shown. The number of time steps required 
for a rotor blade to move by a period was changed from 70 to 
140. The corresponding number of time steps per stator blade 
passing is approximately 133 and 266, respectively. The mean 
pressure showed no change within plotting accuracy. The 
variation of the 1 st and 2 nd harmonic unsteady pressure 
components along the chord is shown in Fig. 3. The unsteady 
pressure magnitude is scaled by the inlet total pressure. The 
nomenclature used is the same as in Ref. [6]: x/c<0 
corresponds to the pressure surface, the leading edge is at 
x/c=0, and x/c>0 corresponds to the suction surface. The 1 st 
harmonic results show very small changes as a result of the 
doubling of the number of time steps per period. The 2 nd 


harmonic pressure plots show small noticeable differences 
between the corresponding results. Reducing the time step in 
half is seen to have only a limited effect on the blade surface 
pressure distribution, while the required computational effort is 
doubled. Based on these comparisons, all remaining 
computations were performed using 70 time steps per rotor 
period. 

The effect of other numerical inputs was also studied. 
Newton sub-iterations are performed at each time step to drive 
the residuals below a prescribed limit. The maximum number 
of such Newton sub-iterations was changed from 3 to 10. The 
changes in the unsteady pressure distributions were smaller 
than those seen in Fig. 3 when the time step was changed. In 
addition, the number of symmetric Gauss-Seidel passes 
performed for each sub-iteration was changed from 3 to 6. The 
changes in the unsteady pressure distributions were extremely 
small as compared to those seen in Fig. 3. Due to the extremely 
limited effect of these changes, the corresponding unsteady 
pressure distributions are not presented. Based on these 
comparisons, all remaining computations were performed using 
3 Newton sub-iterations and 3 symmetric Gauss-Seidel passes. 

Unsteady Flowfields 

The unsteady flowfield for the Run 1 conditions at 85% 
span is shown in Fig. 4 in the form of instantaneous Mach 
number contours at 10 time instants during a passing period. 
The stator wake is clearly seen sweeping across the rotor 
blades. The corresponding flowfields at the 50% span location 
are similar and are not presented. 

Fig. 5 shows a close-up view of the instantaneous 
flowfields at 85% and 50% span. The stator wakes at the 85% 
span location oscillate from side to side, and have a wavy 
appearance. The stator wakes at the 50% span location are 
more uniform in appearance and do not show the unsteadiness 
in location or extent. This spanwise variation in the 
characteristics of the stator wakes is related to differences in 
loading. Further work is required to clarify the change in 
characteristics of the stator wakes along the span. The rotor 
blades at the 50% location show distinct vortex shedding. As 
mentioned earlier, this vortex shedding results in a high 
frequency variation in the unsteady pressure on the rotor blade 
surface near the trailing edge (see Fig. 2c and 2d). In addition, 
the difference in the pattern of shedding between the 85% and 
50% span locations results in a difference in the periodicity of 
the unsteady pressure near the rotor trailing edge (compare Fig. 
2c and 2d). 

Surface Pressure Distributions 

The surface pressure distributions for the Run 1 (small 
gap) and Run 11 (large gap) conditions at 85% and 50% span 
are shown in Fig. 6 and 7, respectively. The experimental data, 
including the error in both variables is also shown for 
comparison. The computational results were obtained using 
two different grids. For the baseline grid, which was generated 
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using established best practices, computations were run on a 
Cray C90 computer due to the large memory requirements; the 
results are labeled as “grid2” in Fig. 6 and 7. A coarser grid 
was used to allow computations to be done on a typical 
workstation used in industry. The results from this coarser grid 
are labeled as “gridl” in Fig. 6 and 7. 

The results calculated on the two different grids, are seen 
to be substantially the same. Some differences are noticeable in 
the phase, but these are typically seen to occur in the regions 
where the unsteady pressure magnitude is small. The CPU time 
required to obtain a converged periodic solution starting from 
an assumed uniform flow was approximately 100 hours on a 
Cray C90 for the baseline grid. The corresponding CPU time 
required was nearly twice as much for the coarser grid on a SGI 
Octane workstation (300 MHz). 

The mean surface pressure is seen to agree quite well with 
the experimental data. In Fig. 6a, the data point near the x/c= 
-0.2 (on the pressure surface) shows a significantly higher 
mean pressure with a very large error band. The computational 
results do not show anything unusual in this location. In Fig. 
6b, the data point near the mid-chord on the pressure surface 
shows a significant drop in pressure, which is also not seen in 
the computational results. It is not clear if there is a 
measurement error or if the computations have failed to capture 
some physical phenomenon. 

The magnitude of the 1 st harmonic surface pressure agrees 
quite well with the measured data. The computational results 
show oscillatory variations along the chord. From the 
measurements at limited number of points, it is not possible to 
say whether or not the variations in the computational results 
are physical. However, results of computations done by the 
authors using the UNSFLO code (not presented here) show the 
same oscillatory nature of the 1 st harmonic pressure magnitude 
and phase. Similarly, the computational results presented by 
Kielb et al. [1] using the UNSFLO code for other conditions 
show the same oscillatory variations along the chord. The 
UNSFLO code is a Quasi-3D Reynolds- Averaged Navier- 
Stokes code. 

In comparing the phase from the TURBO code with data, 
it is noted that the circumferential position of the rotor with 
respect to the stator was not recorded in the experiment. Thus, 
there is an unknown reference in the phase of unsteady 
pressure measurements, which prevents the comparison of the 
absolute phase. However, the variations along the chord can 
still be compared within a constant shift in the phase at ah 
points. The phase from the computational results shows 
reasonable variations except for the (expected) jumps that 
occur when the phase changes sign near ±180 deg. In addition, 
it is noted that the phase presented here is for the adjacent blade 
surfaces of a single passage, rather than for the two surfaces of 
a single blade. 

A comparison of the mean pressure distributions at 85% 
span and 50% span shows that the loading is higher near the tip 
than at mid-span. This is clearly seen by comparing Fig. 6a and 


6b, and Fig. 7a and 7b. This difference in loading level is seen 
in both the data and in the computational results. For the 1 st 
harmonic pressure magnitude, the large peak near the leading 
edge (suction surface) has a higher magnitude at the 50% span 
location as compared to the 85% span location. Similarly, the 
other peaks also have a higher magnitude at the 50% span 
location. These differences, which result from the spanwise 
variation in the stator blade loading, are seen for both Run 1 
and Run 11; see Fig. 6c and 6d, and Fig. 7c and 7d. The phase 
at the two span locations is substantially similar in trends for 
Run 1 (Fig. 6e and 6f) with some differences in the values at 
corresponding chordwise locations, and the chordwise 
locations of some features in the variations. For Run 11 (Fig. 7e 
and 7f), additional differences are seen in the phase. For 
example, the variation seen in phase between x/c=0.0 and 0.2 
in Fig. 7e has no counterpart in Fig. 7f. 

A comparison of the 1 st harmonic pressure for Run 1 
(small gap) and Run 11 (large gap) shows the effect of 
changing the axial gap. Comparing Fig. 6c with 7c, or Fig. 6d 
with 7d shows that increasing the axial gap leads to a reduction 
in the magnitude of the main pressure peak near the leading 
edge. On the other hand, some of the smaller peaks are seen to 
increase in magnitude and the number of the smaller peaks is 
seen to change. The suction surface distributions in Fig. 6d and 
7d illustrate these observations. The effect of such changes in 
the pressure distributions on the overall response is of 
particular interest. 

The results presented thus far have been restricted to the 
Run 1 and 11 conditions. It is neither possible nor useful to 
present ah the results for the other conditions listed in Table 1 
for which computations have been done. Results are presented 
in Fig. 8 showing the mean and 1 st harmonic magnitude of 
pressure at 50% span for Run 4 and 7. These results show the 
effect of the change in rotational speed (19500 rpm to 10100 
rpm), pressure ratio (2.1 to 1.6), and inlet conditions. Also, 
these results show the same good agreement with experimental 
data as seen in previous results. The mean pressure level is seen 
to be higher for the Run 7 conditions and the effect of the 
increased incidence is clear. The 1 st harmonic pressure 
magnitude is seen to have significantly smaller leading edge 
peak for the Run 7 conditions. In addition, the range of 
variation in pressure for the remaining peaks is smaller for Run 
7. 

CONCLUDING REMARKS 

Computations have been performed for a high pressure 
turbine using the TURBO code. The unsteady flowfield for two 
configurations and six conditions has been computed. The 
accuracy of the computational results has been checked by 
varying the relevant numerical parameters. The choice of 
various numerical parameters and inputs is shown to have a 
very limited effect on the results. Results obtained from 
computations on two grids have shown no significant 
differences in unsteady surface pressure. 
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The computational results have agreed quite well with the 
experimental data at both spanwise locations and at different 
run conditions. The mean loading near the tip is seen to be 
higher than at the mid-span location. The trend is reversed for 
the unsteady loading, with the higher unsteady loads occurring 
at the mid-span. A comparison of the surface pressure for Run 
1 (small gap) and Run 1 1 (large gap) conditions has shown that 
the main pressure peak near the leading edge is smaller for the 
larger axial gap. The effect of change in the rotational speed 
and pressure ratio has been shown through a comparison of the 
surface pressure for the Run 4 and 7 conditions. The effect of 
the increased incidence (Run 7) is clearly seen in the mean 
pressure distribution. The 1 st harmonic magnitude of unsteady 
pressure shows a reduction for Run 4 conditions, along with a 
decreased range of oscillations in the values. These effects are 
directly related to the decreased rotational speed. 

Additional work is required to correlate the features in the 
pressure distribution to physical phenomena occurring in the 
flowfield; a space-time plot may be used as a first step in this 
effort. Such work will allow for an understanding of the 
oscillatory chordwise distribution seen in the computational 
results. Future work is also expected to include a comparison 
between the computed and measured blade vibratory strains. 
Additional information regarding the experimental data is 
required before such a direct comparison can be made. In 
addition, future work is expected to include a direct comparison 
of the present computational results with results from other 
CFD codes such as UNSFLO. Such comparison between 
computational results is necessary to verify the results since the 
measurements are available at a small number of pressure 
sensor locations. Further, such a comparison may help identify 
physical phenomena that are captured well by one CFD code, 
or the other, or neither. 
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Table 1: Summary of conditions used in TURBO calculations 


Run 

Mode 

Rotational 

Speed 

(RPM) 

Inlet 

Pressure 

(Pa) 

Inlet 
Temper- 
ature (K) 

Pressure 

Ratio 

Axial 

Spacing 

1 

3 (axial) 

15000 

289590 

296 

2.1 

32% 

11 

3 (axial) 

15000 

289590 

296 

2.0 

46% 

4 

4 (2 nd bending) 

19500 

558495 

553 

2.1 

32% 

16 

4 (2 nd bending) 

19500 

565390 

553 

2.1 

46% 

7 

2 (1 st torsion) 

10100 

289590 

296 

1.6 

32% 

14 

2 (1 st torsion) 

10100 

289590 

296 

1.6 

46% 
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Figure la: Smoothed H-grids used in TURBO calculations 



Figure lb: Close-up view of grid showing details near trailing and leading edges 
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(c) 85% span, near trailing edge 
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Figure 2: Periodicity of unsteady blade surface pressure, Run 1 conditions 
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Figure 3: Effect of number of time steps per period, Run 1 conditions 
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Figure 3 (continued): Effect of number of time steps per period, Run 1 conditions 
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Figure 4: Instantaneous Mach number contours at the 85% spanwise location, Run 1 conditions 
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Figure 5: Comparison of the flowfield upstream of the rotor at the 85% and 50% span, Run 1 conditions 
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Figure 6: Surface pressure distributions for Run 1 conditions 
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Figure 7: Surface pressure distributions for Run 1 1 conditions 
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(a) Run 4 conditions, mean 


(b) Run 7 conditions, mean 




(c) Run 4 conditions, 1 st harmonic magnitude (d) Run 7 conditions, 1 st harmonic magnitude 

Figure 8: Surface pressure distributions for Run 4 and 7 conditions at 50% span 
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